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Abstract 

The goal of this paper is to prove a stable determination of the 
coefficients for the time-harmonic Maxwell equations, in a Lipschitz 
domain, by boundary measurements. 

Introduction 

Let Q be a bounded Lipschitz domain in the three-dimensional euclidean 
space. Assume the medium, modeled by Q, to be non-homogeneous and 
isotropic. Suppose the electromagnetic properties of Q to be described by 
the electric permittivity e, the magnetic permeability \x and the electric con- 
ductivity a. Let E, H denote the electric and magnetic fields, respectively. 
The time-harmonic Maxwell equations at frequency u read 

VxH + iueE = aE, V x E - iufiH = 0, 

whenever the total electric current density is given by aE. Writing 7 = 
e + ia/u, the time-harmonic Maxwell equations can be expressed as 

V x H + iwyE = m 
VxE-iu/iH = 0. ^ ' 

It is known that this system, complemented with a suitable prescribed data 
on the boundary, is well-posed for some w's. In fact, we have the following 
result. 
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Theorem 1 Let Q be a bounded Lipschitz domain and fi,e,a G L°°(Q) 
satisfying 

//>//> 0, e > e' > 0, tr > 0; 

a. e. in f2, with //',£' positive constants. Given T G TH(dQ), the problem of 
finding E,H G iJ(fi;curl) solving ([T]) in and satisfying either N x E = T 
or NxH = T is well-posed for any wgC \ {0} except for a subset of 

{u G C : - ||cr/£|| L oo( Q ) < Imw < 0} 

with no accumulation point in C \ {0}. 

Here N stands for the unit vector field normal to dQ, the boundary of Q. 
Precise definitions of the spaces H (Q; curl) and TH(dQ) are given in Sub- 
section 11.11 The frequencies u for which the direct problem is not well-posed 
are called resonant frequencies. 

A proof of this result in the case where dQ is C 2 can be found in [29J. 
The well-posedness of the same problem in the context of non-smooth do- 
mains is studied in [H] but the result stated there assumes a to be zero. 
Finally, a proof of the precise statement of Theorem [T] can be given using the 
compactness result stated in [33] and a lemma due to Peetre (see [23J and 

my 

Theorem [T] allows us to model the boundary measurements by means of 
the admittance or the impedance maps. The admittance map is defined as 

A ad : T G TH(dtt) i — ► NxH G TH(dtt), 

where E, H is the solution for ([1]) with N x E = T. The impedance map is 
defined as 

A im : T G TH{dVt) i — y NxE G TH{dVt), 

where E, H is the solution for ([1]) with N x H = T. It is a consequence of 
Theorem [T] that these maps can only be used out of resonant frequencies. 

With one of these maps at hand, one can set the problem of determin- 
ing the coefficients of equations in ([T]). This inverse problem was initially 
proposed by Somersalo et al in [2H]. In these notes we prove not only deter- 
mination of the coefficients but also stable determination. 

When trying to determine in a stable manner the coefficients by one of 
these maps, we have to face the problem of choosing u > to be non-resonant 
for the class of coefficients to determine. How can we manage to solve this if 
the position of the resonant frequencies depends on the unknown coefficients? 
As it was pointed out in [21], the same happens in practice. How can we 
know when our choice of frequency is close to a resonant frequency? 
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In order to avoid this problem we shall model the boundary measurements 
by the Cauchy data set. Given a frequency u > 0, the Cauchy data set is 
defined as follows: (T, S) G C(/i, 7 ) if and only if (T, S) G (TH(dn)) 2 and 
there exists a pair (E, H) G {H(Q; curl)) 2 solution of ([1]) satisfying NxE = T 
and NxH = S. 

This way of modeling the boundary measurements has been used success- 
fully in [5], [21] and [22] • However, as far as the author knows Cauchy data 
sets have not yet been used to establish any result of stability. 

In order to measure the proximity of the Cauchy data sets associated to 
given pairs //i,7i and /X2,72 we introduce a pseudo-metric distance. 

Definition 1 Let /xi,7i and fi 2 ,j2 be two pairs of coefficients. Consider u 
a positive frequency and let Cj denote Cipj^j). Let us define the pseudo- 
metric distance between the Cauchy data sets C\ and C 2 as 

&c(C u C 2 ) = max sup inf || (7), 5j) - (T k , S k ) \\ (THm))2 ■ 

ll 7 fcllTH"(an)) =1 

We say that 5c is a pseudo-metric distance because if 5c{C\, C 2 ) = 0, we can 
only ensure that C\ —C 2 . 

The definition of 5c is inspired in the Hausdorff distance. Unlike the 
latter distance the former one is comparable to ||A" d — A| d I when uj is a non- 
resonant frequency for /ij,7j with j = 1,2. Here |.| denotes the operator 
norm for linear operators on TH(dQ). 

Our result requires certain stability of the problem on the boundary and 
since this has not been proven yet, we shall introduce some definitions. 

Definition 2 Given two constants M, s such that < M, < s < 1/2, 
we shall say that the pair of coefficients /i, 7 is admissible if they satisfy the 
following conditions. 

(i) Uniform ellipticity condition. The coefficients 7,/i G C 1,1 (fi) satisfy 

M~ l < Re7(x) M _1 < fi{x); 

for any iGSl. 

(ii) A priori bound on the boundary. The following a priori bound holds on 
the boundary 
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(iii) A priori bound in the interior. The following a priori bounds hold in 
the interior 

Il7llw 2 .°°(tt) + ll^llw 2 .°°(«) — ^> II7|Ih 2 + s (q) + ll^llffa+«(n) — M. 

Definition 3 Let M, s be the constants given in Definition [2] and let u be 
a positive frequency. We shall say that a pair /i, 7 is in the class of B -stable 
coefficients on the boundary at frequency u if /x, 7 is an admissible pair and 
there exists a modulus of continuity B such that, for any other admissible 
pair ft,, 7, one has 

117 - 7llcfo.i(«i) + II/ 1 - ^Wco.naa) ^ B ( 6 c(C, Cyj , 

l|V(7-7)llx~(«i;C3) + HV^-^Hioc^s) < B (6 C (C,C)) . 

Here C, C are the Cauchy data sets associated to the pairs /x, 7 and jl, 7, 
respectively. 

With these definitions at hand the stable determination of the electro- 
magnetic coefficients can be stated as follows. 

Theorem 2 Let Q be a bounded Lipschitz domain and let u be a positive 
frequency. Consider /Ji,7i and /i2,72 any two pairs in the class of _B-stable 
coefficients on the boundary at frequency uj, with B satisfying \r\ < B(\r\) for 
all |r| < 1. Then, there exists a constant C = C(M) such that the following 
estimate holds 

||7i -72|lHi ( n) + llA*i — A*a II ffi(n) ^ C\\ogB{5 c {C 1 , C 2 ))|" A , 

for some constant A such that < A < 2/3s. Here Ci,C 2 are the Cauchy 
data sets associated to the pairs jUi,7i and /i2,72, respectively. 

As in the inverse conductivity problem, it should be possible to prove that any 
admissible pair is in the class of Holder-stable coefficients on the boundary 
for any frequency u, that is, with -B(|r|) = \r\ a for < a < 1. Notice that 
in the conductivity case a logarithmic module of continuity, as the one in 
Theorem El is optimal (see [16J). 

In [20] and [22] the coefficients were assumed to be constant on the bound- 
ary -which is actually quite natural from the point of view of applications-, 
in that case our result reads as follows. 

Corollary 1 Let Q be a bounded Lipschitz domain and let u be a positive 
frequency. Consider /ii,7i and /i2,72 any two pairs of admissible coefficients 
-in the sense of Definition [2j Assume that 

fallen = A*2|an, d x it><i\dn = d^^ldn, 71 1 an = l2\dn, d x jji\QQ = d^^lan, 
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with j = 1,2,3. Then, there exists a constant C = C(M) such that the 
following estimate holds 

Il7i-72|lifi(n) + 11^1-/^11^1(0) ^ C\\og5 c (C 1 ,C 2 )\' x , 

for some constant A such that < A < 2/3s. 

When the boundary is smooth, at least C 1,1 , it is possible to prove the 
existence of H ^solutions, this allows us to describe traces not only in a weak 
sense but also in a strong sense. In consequence, the boundary measurements 
can be modeled in better spaces (see [20] and [22]). Furthermore, it is well 
known (see [26], [27] [3], [I] and [8]) that Maxwell's equations may not admit 
i7 1 -solutions even with if 1//2 -boundary data, whenever the domain is neither 
convex nor has C 1 ' 1 -boundary When working on Lipschitz domain, this 
lack of smoothness makes necessary to introduce some non-standard Sobolev 
spaces. 

The unique recovery of C 3 -coefficients 7 and fi from boundary data was 
proved in [20], and later simplified in [22], in the context of C 1 ' 1 -domains. 
The regularity in our result of stability agrees with the required in [20J and 
[22], nevertheless from our proof one can state a uniqueness result for C 1,1 - 
coefficients in the context of Lipschitz domains. Boundary determination 
results were given in [T7J and [12] in the case that the boundary is smooth. 
The more general chiral media was studied in [TS] . For a slightly more general 
approach and more background information, see also the review article [2"T| . 

In [28] the relation between the inverse scattering problem (for short 
ISP) and the inverse boundary value problem (for short IBVP) was studied. 
Further, Sarkola proved that the former problem can be reduced to the latter 
one. On the other hand, in [9 J Hahner gives a proof of the stability for the ISP 
when magnetic permeability /i is constant and as a consequence he obtains 
stability for the IBVP for constant magnetic permeability when the domain 
is a ball. Thus, our corollary generalizes the stability result stated in [9], at 
least, for the IVBP. 

The main motivation to study the inverse problem in the setting of Lip- 
schitz domain comes from the technical hypothesis that had to be assumed 
in [0] when studying local data result. Some of the ideas shown along these 
notes can be adapted to prove a stable determination of the coefficients from 
local data. This is accomplished in [7J. 

The structure of the paper is as follows. There are two sections and an 
appendix, in the first section some preliminary details are given, while in 
the second one we prove the stability of the IBVP. The appendix contains 
basic definitions and properties of the spaces used along these notes. In a 
deeper extent, the first subsection of Section [1] is dedicated to the functional 
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spaces adapted to Maxwell's equation and their traces. The definitions of 
these spaces and traces have not been placed in the appendix in order for 
the reader to get used to them. In the second subsection of Section (TJ we 
transform the system (TTJ) into a Schrodinger-type equation. In Section |2] we 
perform some of the standard steps when studying this IB VP: prove a suitable 
estimate relating the boundary measurements with the coefficients in the 
interior, construct special solutions for the Maxwell's equations through the 
Schrodinger-type equation and plug these solutions into the suitable estimate. 
Finally, we obtain the estimate stated in Theorem [2] by the use of a Carleman 
estimate. 

Acknowledgement. This paper is part of the author's doctoral disserta- 
tion and it has been written under the supervision of Alberto Ruiz. The 
author would like to thank him for his support -gracias por dejarme elegir y 
corregir mis errores. Additionally, the author would like to thank the com- 
ments of the referees. The presentation of the result has improved a lot with 
their suggestions. The author was economically supported by Ministerio de 
Ciencia e Innovacion de Espana, MTM2008-02568-C02-01. 

1 Preliminaries 

Along this section we introduce some spaces that turn to be useful in the 
context of Maxwell's equations and we transform ([1]) into a Schrodinger-type 
equation. 

Let us denote by E the three-dimensional euclidean space and by TE the 
module of smooth vectors fields over the real smooth functions. The elements 
u G TIE will be called real vector fields and can be expressed as 

U = ( u (2) u (3) )*. 

The euclidean metric induces on E a volume element dV and the euclidean 
distance d e . Additionally, for any two vector fields u, v, we shall denote by 
u ■ v and u x v the point-wise inner product and the point-wise cross product 
of vector fields, respectively. On the other hand, denote 

XE = {u + iv : u,v G TE}. 

The elements of XE will be called complex vector fields. 

Finally, the gradient, divergence and curl operators will be denoted by 
V, V- and Vx, respectively. 
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1.1 Non-standard Sobolev and Besov spaces 

Most of the facts collected here can be found in j 1!) . 

Definition 4 Define the spaces 

H(tt; div) = {ue L 2 (VL] C 3 ) : V • u G L 2 (tt)}, 
H (fi; curl) = {t)6 L 2 (ft;C 3 ) : V x t> G L 2 (fi;C 3 )} 

equipped with the graph norms 

IMIiT(fi;div) = H M llL 2 (fi;C 3 ) + ' U Wl 2 (Q.) ' 

IMItf(fi;curl) = IML 2 (fi;C 3 ) + \\^ XV \\L 2 (Q;C :i ) ■ 

Now we define the traces of elements belonging to these spaces. For any 
u G H(Q;div) the normal trace of u, that is, A" • u, can be defined as an 
element of B~ l l 2 (dVl). Namely, for any g G B 1 / 2 (dVt) 

(N-u\9)= [ iy-u)JdV+ [ u-VfalV, (2) 
Jn Jn 

where / G H l {Vl) and f\dti = g- It is well-known that 

N ■ . : H(Q; div) — > B- 1/2 (dQ) 

is a bounded and onto operator. Additionally, note that if u G H(Q;div), 
f G CP'^dQ) and f is any extension of / such that / G C - 1 ^) (see §U\ for 
extensions from closed set), then 

(fN-u\g) = (N-u\gJ) = (N-{fu)\g) , (3) 

for any g G B 1 ' 2 (dVL). 

Similarly, for any u G H (fi; curl) the tangential trace of u, that is, Nx 
u, can be defined as an element of B~ l l 2 (dVL; C 3 ). Namely, for any w G 
B l l 2 (dVL] C 3 ) 



(A^xn|u;)= (Vxu)-vdV - / u-(Vxv)dV, 
in in 

where t> G C 3 ) and u|gn = w. In this context, the operator 

Nx.: H(Q;cm\) — > B~ 1 ' 2 (<9ft; C 3 ) 
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is bounded but not onto. Let its range be denoted by 

TH(dn) = {w G B- l/2 {dVt; C 3 ) : 3u G H(Q; curl), Nxu = w}. 
This vector space can be equipped with the norm 

\tH(90) = inf {ll M L(f7;curl) : U E #(ft;CUrl), N X U = w} , 



\w\ 



which makes it a reflexive Banach space. Moreover, TH(dQ) is continuously 
embedded into B~ 1 ^ 2 (dQ; C 3 ). In addition, the map 

Nx.: TH(dtt) — > (TH(dtt))* 

which, for w±,W2 G TH(dfl), is given by 

(Nxw l \w 2 }= / (Vx«)-udV- / u-(Vxt)) dV, 
Jn Jn 

where u, v G curl) are such that iVx-u = u»i, Nxv = W2, is well-defined, 
bounded and an isomorphism. In particular, one has that (iVx.) -1 = —Nx., 
(Nx.)* = -Nx. and 



({Vxu)-vdV = [ u-{Vxv)dV - (Nxu\Nx(Nxv)) . 
Jn Jn 



(4) 



Before going further let us point out that, ifu> G TH(dtt), f G C°' l (dQ) 
and / is any extension of / such that / G C 0,1 (f2), then 

(fw\z) = (w\Jz) = (Nx(fu) 

for any z G B 1/2 (dQ; C 3 ) and any u G H(Q; curl) such that Nxu = w. Note 
that last pairing implies that fw G TH(dQ). 

Another identity to mention is the following, if w G TH(dQ), f G 
C 0,1 ((9fi) and / is any extension of / such that / G C 0,1 (fl), then 

(fw\Nxz) = (^Nx(fu) Nxz^ = (w\Nx(fz)), (5) 

for any z G TH(dQ) and any u G if (Q; curl) such that N xu = w. 

We can define the surface divergence operator over the space TH(dQ) 

Div : TH{dtt) — > B- 1/2 {dQ), 



as 



Divw = -JV-(Vxu), 
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where u G H(Q; curl) and Nxu = w. Since Vxm G H(Q; div), -IV-(Vx?i) 
makes sense and it belongs to B^^^dVt); this operator is well-defined and 
bounded. 

The surface divergence can be used to give an intrinsic description of the 
space TH(dQ). In [19] Mitrea proved that there exist constants C\,Ci > 
such that, for any w G TH(dQ), the following estimates hold 



\w 



\ W \\TH(dn) ^ ^2 



MI.B-i/2(af};c 3 ) + ll Divu, llB-i/2(an) / • 



(6) 
(7) 



Lemma 5 There exists a positive constant C such that the following esti- 
mate holds 

< C 



\w\ 



TH{dQ) 



\\f W \\TH(dn) — ° \\J HC ' 1 ^) 

for any w G TH(dn) and any / G C ' 1 ^). 

Proof: Consider / an extension of / such that / G C°' l (Q) satisfying 



/ 



< C 



and any u G if curl) such that N-u = w. Then 



ll/H 



TH(dn) 



< 



fu 



< 



H(n;curl) 



'U 



_ff(n ; curl) 



< c 



\C°> 1 (dU) ll M ll_ff(C;curl) 

Taking infimum in u G curl) one gets the estimate. 



□ 



1.2 Maxwell's system as a Schodinger equation 

In this section we shall transform Maxwell's equations into a Schrodinger- 
type equation. The idea of this transformation was already introduced in 

[22]. 

There is a well known process which allows us to transform Maxwell's 
equations into a Schrodinger-type equation, in order to do so we require some 
extra smoothness of the coefficients, namely /x, 7 G C 1,l (VL). The first step in 
this process is to augment the Maxwell system with two scalar equations: 

V ■ (jE) = 0, V • (jjtH) = 0. 

The information coded in these scalar equations was already present in the 
initial system. In order to check this, it is enough to take divergence in each 
equation in ([!]). 
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Next, we introduce a new system inspired in the four equations mentioned. 
This new system reads as 

V • (7.E) + iwyfih = 
-7 _1 V(7e) + V x E - iou/iH = 

V • (fJ>H) + iujfxe = 
/j _1 V(/j/i) + V x H + iu-fE = 0. 

The new terms preserve the physical units of measure of the original four 
equations. The new system, called henceforth augmented system, can be 
written as matrices in the following way 







D- 


\ 




D 


-Dx 




D- 








D Dx 






/ 



+ 



( Ufl 




Da- 






Da 




Dfi- 




V m 




wyl 3 



( h \ 

H_ 

e 

\ E ) 



0, 



where a = log 7, j3 = log/x, Ij is (j x j)-identity matrix, with j e N and 

1 



D- = - ( d x i d x 2 d x 3 ) , 



D 





|. 




\ d xi J 





d x 3 
—d r 2 d T i 



-d x 3 d x 2 
-d x i 



In a much more compact manner we shall express the augmented system as 
(P + V)X = 0, where 



P = 







D- 


\ 




D 


-Dx 




D- 








D Dx 






I 



( UjJL 



V = 



DP- 



\ DP 



Da- \ 



Da 



W7/3 / 



Note that E, H is a solution for Maxwell's equations, if and only if, X 1 = 
( h H f I e E t ) is a solution for the augmented system and the scalar fields 
e, h vanish. 

The next step is to rescale the augmented system, that is 



(P + V) 
where 



ti- 1/2 h 






T ir2 h 



Y 



T l/2 h 









(P + W)Y, 



W = K Is + 





/ 




Da- \ 


1 






Da Dax 


2 




DP- 






V DP 


-DPx 


I 



(9) 
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with k = o;/i 1//2 7 1//2 . We shall call 

(P + W)Y = 

the rescaled system. 

The advantage of rescaling is that 

= (P + W){P - W t )Z = (-A/ 8 + Q)Z, 
= (P - W l )(P + W)Z' = (-AJ 8 + Q')Z' 
= (P + iy*)(P - W)Z = (-A/ 8 + Q)Z, 



(10) 

(11) 

(12) 



where Q,Q',Q are zeroth-order terms. Here W t denotes the transposed of 
W and W* stands for W t . No first order terms appear in ffTOl , ( TTTj) and f[T2"j) , 
giving as a result a Schrodinger-type equation. Mind 

Q = -PW l + WP- WW 1 . (13) 

Note that if Z is a solution for ffTO]) in f2, then F = (P — W l )Z is a solution 
for the rescaled system in Q, hence 

/i" 1/2 / 4 



X 



1 



is a solution for the augmented system. In the same manner, if Z is a solution 
for (H2D, then F = (P - is a solution for (P + W*)y = in ft. 
For later uses, 





/ Aa 


\ 


1 


2V 2 a - AaJ 3 


2 




A/3 




V 


2V 2 /3- A/3/ 3 / 



(k 2 + l(Pa-Pa))/ 4 



2Pk- 



2P/t 



2Dk- 



2Dk 



(n 2 + \(D(3-D(3))h 



(14) 



/ 



Q' = ~: 





/ A/3 


\ 


1 


2V 2 /3 - A/3/3 




2 




Aa 




V 


2V 2 a-Aa/ 3 / 



/ 
V 



(k 2 + 1(P/3-P/3))/ 4 


\ 

2Dkx 


-2Dkx 


(k 2 + \{Da-Da))h 



(15) 
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and 





/ -A/3 


\ 


1 


-2V 2 /3 + A/3J 3 




2 




-Aa 




V 


-2V 2 a + AaJ 3 / 



(k 2 -1(DP-DP))I 4 


\ 

-2D~Kx 


2Dkx 


(k 2 - \{Da-Da))h ^ 



(16) 



Q 



/ 
\ 

with V 2 / = (C, /);.„ ,• 

The computations needed to get ()14p and ffT5l) can be found in [6]. The 
same kind of computations gives (fT6l) . 

In order to make as concise as possible the presentation of our proof, we 
introduce some additional notation. Let Y, Z be in the form 

Y={f 1 (u 1 )* 1 f 2 (u 2 yy, Z={g' (viy\g 2 {v 2 f )* , 



define 



2 

(Y\Z) n = [Jj J 3 JdV + j/^ dV 

(y\z) 9 n = J2([ f j 9 JdA + [ 

• =1 \Jdn Ja 



u J -vidA . 



In the first identity we are assuming f^gi G C°°(Q) and u\ t> J G Af]E|^j with 
j = 1,2, while in the second identity p, g j G C°°(<9f2) and u j ,v j G A'Elan 
with j — 1,2. The following integration by parts holds 

(PY\Z) n = {P N Y\Z\ m ) m + (Y\PZ) n . 

Here, when A is a (possibly complex) vector field we denote 

/ 

F\ - , \ A ~ AX . (17) 



Finally, for elements Y in the form given above we define, for \s\ > 0, 







A- 


\ 




A 


-Ax 




A- 








A Ax 






J 



\Y\ 



E (11/ 



+ \\u J 



l_H" s (n ; c 3 ) 



i=i,2 



12 



and 



\ Y \\L''(n-y) 



E (11/ 



3=1,2 



\L 2 {n) ^11" llL 2 (f7;C 3 ) 7 ' 



Here and henceforth, ^ denotes the space C x C 3 x C x C 3 . On the other 
hand, we define, for < \s\ < 1, 



and 



\Y\ 



\Y\ 



B°{dn-,y) 



( ll^ lls s (an) ll^ ll ;; ,,«>- ! i / ' 



ls s (an ; c 3 ) 



3=1,2 



E (11/ 



L-aiv.-.y) - y\\J 3 \\L^(an) + IMIl^ck^c 3 ) 

3=1,2 



2 A log-type estimate 

In this section we prove Theorem [21 The proof consists of three ingredients, 
namely, a suitable estimate relating the electromagnetic coefficients in Q with 
their corresponding Cauchy data sets, the constructions of special solutions 
and the use of a Carleman estimate. 

Lemma 6 Let fij, jj belong to C 0,1 (fi). Then one has that, for any Y\ given 
by 

n = (o ^|o 7^)' 

with Ei, Hi G i/(fi;curl) solution for jT]) in with coefficients /ii,7i, and 
any 

Y 2 = ( f 1 (w 1 )' I / 2 {u 2 Y Y e H\Q) x H(Q; curl) x if 1 ^) x H(Q; curl) 
solution for (P + W 2 )Y 2 = in f2; the following estimate holds: 



+ 



<C5a(Cx,C 2 ) 

1/2 



/'2 



1/2 



72 



c - 1 (an) 



S , i|lB 1 /2(an) + 



I #2 1 1 £i/2(,9 n ) 



1/2 
72 



\ Z l\\TH(da) 



+C7(||JVxS 1 || rH(an) + 11^x^1 



1/2 

M 2 



lT_ff(9Q) 



-1/2 -1/2 



1/2 1/2 



+ 



1/2 1/2 

7i -7 2 



-1/2 -1/2 

7i -7 2 
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Here 5-1,5-2 G P 1/2 (<9fi) stand for g x = f^gn, g 2 = Plan while z x ,z% G 
TH{dQ) stand for % = Nxu 1 , z 2 = Nxu 2 . Here (7,-, with j = 1,2, stands 
for the Cauchy data set corresponding to fij, jj. Recall that W 2 is the matrix 
dSD associated to ^2,72- 

Along these notes, /ij and 7, should be understood either as themselves 
or as their traces, according to the context. 
Proof: Let L be 



L=( ^ 2 HI\0 ^El) 



with E 2 , H 2 G H(Q; curl) an arbitrary solution for (pQ) with coefficients /i 2 ,72- 
Since (P+H/" 2 *)y 2 = and (P+W 2 )£ = 0, one has that (L\PY 2 ) Q = (PL\Y 2 ) Q , 
hence 



(Y 1 \PY 2 ) n - (PYi|y 2 ) n = (Fi - L|PF 2 ) n - (P(F! - L)\Y 2 ) n . 
On the other hand, we have, using (J2J), (j3J), © and (jSJ), that 
(Y 1 -L\PY 2 ) n -(P(Y 1 -L)\Y 2 ) n = 
= i(N-(fi 1 H 1 - fi 2 H 2 ) fi 2 1/2 g 2 ) +i(N-(ji 1 H 1 ) (/i7 1/2 - fi 2 1/2 )g 2 



+1 (iV-(7iPi - 72P2) 7 2 1/2 9i) + * (iV-( 7 iPi) (7r V2 - 7 2 _1/2 )^i 
-i (nx(H! - H 2 )\Nx(ii 1 2 /2 z 2 )^ - % (nxR^Nx^ 2 - fil /2 )z 2 )^ 
+1 (nx{E 1 - E 2 )\nx (72 /2 ^i)) + i (NxE^Nx^l 12 - 7 2 1/2 )^i)) • 
Furthermore, from Maxwell's equations one deduces that 

N-h.EA = — DivfiVxFA iV-(u,F 7 ) = —DivfNxEA, 
ioj iuj 

for j = 1,2. Hence, denoting NxEj = Tj and NxHj = Sj we obtain 



(Y 1 \PY 2 ) n -(PY 1 \Y 2 ) n = 
-- (Div(T 1 - T 2 ) fi 2 1/2 g 2 ) (DivTi (>7 1/2 - fx 2 L,z )g 2 



-l/2^ 



+- (Div(5i - 5 2 

u; \ 



72 _1/ ^i) + -(Div5 1 



-% (S 1 - S 2 
+1 (Ti - T 2 



(7i 1/2 -7 2 1/2 )^i) 
Nx{^ 2 z 2 )) - % (s^Nx ((/ij /2 - /if )* 2 ) 



iVx (7^)) + i (TxjiVx (( 7l 1/2 - 7 2 1/2 )^i)) • 
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By using the appropriate dualities and the estimates 

\(Y 1 \PY 2 ) a -(PY 1 \Y 2 ) a \< 



and (jSJ) we get 



< C (\\T X - T 2 \\ TH{dn) + || - S 2 \\ TH{an) 



1/2 
72 



INI 



TH(dn) 



+ 
+ 



1/2 
M 2 



1/2 1/2 

7i -7 2 



\\ z i\\TH(dn)^ + C 
"zi" 



72 



1/2 



-1/2 
Ml 



-1/2 
M 2 



-1/2 
M 2 



\\9iW 



B 1 / 2 (dU) 



s 1 /2(an) 



+ 



1/2 

Mi 



1/2 
M 2 



TH(dn) n 



7i 



-1/2 



-1/2 
72 



INI 



IITill 



TH(an) + ll^i|lTi?(9n) 



B 1 / 2 (an) 



This estimate holds for all (T 2 , S2) € C 2 , since £"2, i?2 was chosen to be an 
arbitrary solution for (pQ) with coefficients /i 2 , 72- Finally, the wanted estimate 
is a consequence of Definition [TJ □ 

Proposition 7 Let 71, //1 and 72,^2 be in the class P-stable coefficients on 
the boundary at frequency u, with 5 as in Theorem [2J Then, there exists a 
constant C(M) such that, for any Z x G iP^O; 3>) satisfying Yi = (P X -W{)Z X 
with Yi as in Lemma |6] and any Y2 £ H 1 ^; y) as in Lemma El one has 

I («?! - £2)^2)0 I < CB{8 C (C X , C 2 )) ll^iH^^) r 2 ||flx (niy) 

+ C , S(5 C (Ci,C 2 )) (||-&l|| H (n;curl) + \\ H 1 II H(n;curl) 



I/ 1 ! 



Hi(n) 



li?(fi;curl) 



lff(fi;eurl) 



+ II/ 2 



lff!(n) 



Here is the matrix (|T3"|) associated to /ij, 7,- with j = 1,2. 
Proof: From (1131) one has 



((Qi - Q 2 )Zi|y 2 ) n = - (P(w? - Wl)^i|y 2 ) n 

+ ((Wi - W 2 )PZ 1 |y 2 ) n - ((WiWf - WW^Y^ 
= (W - W t 2 )Z 1 \P N Y 2 ) m - (W t 1 Z 1 \PY 2 ) Q + (iy 2 %|PY 2 ) n 

+ (Wi{p - w{)z x \Y 2 ) n - (pz 1 \w;y 2 ) q + (W 2 t Z 1 \W 2 *Y 2 ) Q 

= {(W* - W 2 t )Z 1 \P N Y 2 ) m + ((P - W t 1 )Z 1 \PY 2 ) Q + (W 1 (P - W t 1 )Z 1 \Y 2 ) n 

= {{w[ - w 2 t )z 1 \p N Y 2 ) dn + (yi\PY a ) a - (py a |y a ) n . 

In order to get the penultimate identity, we used twice that (P + W 2 )Y 2 = 0. 
In the last one, we used that Y X = {P- W{)Z X and that (P + W 1 )Y 1 = 0. 
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It is a straight forward computation to check the next estimate 
| {{W{ - W 2 t )Z 1 \P N Y 2 ) m | < C (|K - K 2 \\ Loom + ||V(/?i - h)\\ L °°(an ; c 3 ) 
+ \\V(at - tt 2 )|| LK > ( a n;C 3)) II^iII^^q^) \\y2\\ L 2 {m . y) ■ 
Here, as usually, the norm of L°°(dQ; C 3 ) is 



\w\ 



L°°(an ; c 3 ) 



i=i 



\w 



(i)| 



iL°°(an) 



for any w G <Y(E)|gn. 

It is a routine computation to check that, on one hand 

IK - K 2 \\ LX{m) < C (J|7 2 - 7i|| L <»(an) + Wto ~ ^illi-(sn)) 
<CB(6 C (C U C 2 )), 

°°(dn ; c 3 ) °°(dn;C 3 ) 

<CB(6c(C u C 2 )), 

°°(dn;c 3 ) ^ ^ (ll/^i ~~ ^ILoo^n) + II ^(^1 ~ ^2)|| L oo( 9Q;C 3) 
<CB(8 C {C U C 2 )). 

and on the other hand, 



fh 



1/2 



co. 1 (an) 



+ 



72 



1/2 



co.^an) 



+ 



1/2 

M 2 



co.^an) 



+ 



1/2 
72 



c .!(9n) 



< C 



-1/2 -1/2 
-1/2 -1/2 

7i -7 2 



1/2 1/2 
Ml 

1/2 1/2 

7l -7 2 



co-Han) 



c - 1 (an) 



co. 1 (an) 



<C\\^-^\\ a ^ m <CB{S c {C u C 2 )), 
<C||7i-72|lao, H an)<C^(Wi>C2)), 
<C|| A *i-/i 2 || C o.i (fln) <C7B((y C r(Ci l C7 2 )) l 
<C\\j l -'y a \\ ao , 1(ga) <CB(S c (C u C 2 )). 



c ' 1 (an) 

Putting together all these estimates and Lemma El we get 

I ((Qi - Q 2 )^i|y 2 ) n | < CB{8c{C x , C 2 )) H^U^an^) \\Y 2 \\ Bl/2{dn . y) 
+ CB(5 C (C 1 ,C 2 )) (\\NxE 1 \\ THidU) + \\NxH 1 \\ TH{m) ) 

x (llfl'i|lBV2(an) + ll z illrH(an) + llfl^Hsvatan) + H^Hth^) 
hence we deduce the estimate given in the statement. □ 



16 



2.1 Construction of special solutions 

Here we construct two kinds of special solutions, one for the Schrodinger- 
type equation and another one for (P + W*)Y = 0. The first one was 
already constructed in [22J but we give here the proof in order to keep track 
the constants. The second kind of solution is inspired on the solutions given 
in P3]. 

Let B(0; p) be the open ball centered at the origin O with radius p > 
and such that Q C B(0;p). Sometimes B(0;p) will be denoted by B p 
to simplify the notation. Let 6$ and po denote the electric and magnetic 
constants, respectively. Extend the coefficients 7, p defined in Q to functions 
in E -still denoted by 7, p-, preserving their smoothness and in such a way 
that 7 — e , p — po have compact support in B(0; p) (regarding to extensions 
see [30J). Note two simple facts. Firstly, the extensions still satisfy the 
a priori bound and the a priori ellipticity condition in E. Secondly, the 
extensions of the matrices (JHJ), ({15]) (JIB]) -still denoted by Q,Q',Q- satisfy 
that u 2 e p I 8 + Q, u 2 e Q p I 8 + Q' and u 2 e p Q I 8 + Q have compact support 
mBiO^p). 

We shall construct solutions for ffTOl) in E with the form of a complex 
geometrical optic solution (CGO solution for short), that is, in the form 

Z = e iC - x {L + R), 

with L = L(C) constant and £ G C 3 . 

Lemma 8 Let S be a constant such that —1 < 5 < and let £ e C 3 be such 
that C ' C — l ^ 2£: oA i o with 

\c\>cj2\\(^oh+Q)n L ^ Bp) 
j,k=i 

for some constant C = C{5, p) > 0. Then, there exists a 

Z = e Kx {L + R) 

solution for (-AI 8 + Q)Z = in E, with Z\ a G # 2 (ft; 3>), L = L(() constant 
and R = R(() satisfying 

g 

m w- S IK w W8 + Q)?|| LOO(Bp) , (19) 

8 

||P^||^ < C(^P) E ll( w2e o^8 + G)ilU ( B P ) + ll^ll^y)- ( 2 °) 

j,k=i 
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The norm in the lemma is 



L|(E;C3) ; ' 



i=l,2 



with y = ( f 1 (u 1 ) 1 | / 2 (w 2 )' ) and 

\\f\\ll = [ (l + \x\ 2 ) x \f\ 2 dx, 

for < |A| < 1. 

Proof: It is an easy matter to check that 

e-^ x (-Ah)e iC - x (L + R) = [-A - 2i( ■ V + C • (]h(L + R), 

hence, if Z = e l ^' x {L + R) is a solution for (IIP]) , then i? solves 

((-A - 2< • V + u 2 e 0f x )h + Q)i2 = -(u 2 e 0f i h + Q)L. (21) 

We will use this equation as the starting point for the construction of the 
CGO solutions. 

It is a very well known fact that the operator -standing for the convo- 
lution with the fundamental solution for Fadeev's operator (—A — 2i( ■ V)- 
with | C| > 1 satisfy the estimate 

IKVM^f 11/11^, (22) 

Rjcyiu < c{6) \\f\y , (23) 

for any / G £ 2 +1 with j = 1,2,3. Estimate fT22|) was proven in [31], while 
fT23|) was proven in [2]. 

Denote = G^I 8 - Applying to both sides of f[2"Tj) we get 

(J 8 + F c (u 2 e fi I 8 + Q))i? = -F ? ( W 2 ^o/8 + Q)L. (24) 

On the other hand, we can estimate 

\\F ( (uj 2 e ^l8 + Q)R\\ L p < S \\(" 2 £oVoh + Q)j\\ L oo (Bp) \\R\\i? s y, 

j,k=i 

where we applied (J2"2j) and used the fact that u 2 e [i Ig + Q has compact 
support in B(0;p). 
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Then, assuming 

8 

^ rifx n \ ST^ \\f, ,2„ ,, j i n\k\\ 

\L°°(B P ) 



|C|>C(ME \\{u 2 e ^h + Q))\ 
j,k=i 

the operator (J 8 + F^(cu 2 £oA i o-^8 + is bounded in L 2 y and 

i? = -(/ 8 + F c {uj 2 e ^h + Q)Y 1 F c {u 2 e Q ^h + Q)L, 

with 

g 

W R hp < [7H l L l Z ||(cAo/^8 + Q)i|| Loo(Bp) • 
M i,k=l 

Here we used again fl22l) and the fact that oj 2 £:oA i o^8 + <5 has compact support 
in B(0; p). This compactness is crucial in our arguments. 
On the other hand, from (EMI) and (|2"5|) we deduce that 

II^ILjy < ||P[F c (a; 2 £o/io/8 + Q)(i2 + iv)]|| L , 3; 

8 

j,k=i 

□ 

The arguments given in the above proof were used in |31j for scalar equa- 
tions of the same type as (12T1) . We exploit below the uniqueness for these 
scalar equations. 

We will use the CGO solutions constructed for (TTDT) to produce analogous 
solutions for Maxwell's equations. The procedure follows the ideas exposed 
in Section [L~2| using the decoupled scalar equations in ffTTl) . 

Proposition 9 Let 5 be a constant such that —1 < 5 < and let £ e C 3 be 
such that C " C — k^o/^o with 

|C|>C(M \J2 ll w2e °**o + gjL« (B) + Z \\(u 2 eojJ h + Q) k 



L°°(B P ) 1 II V~ -uA-u-o i -Wj || L °°( Bp ) 

J=l,2 j,k=l 



where 



qi = - k 2 - - A {Df3-Df3), q 2 = - k 2 - -^Da-Da). 

Then, there exists a 

Z = e^ x (L + R) 
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solution of (-AI 8 + Q)Z = in E, with Z\ a G H 2 (Q; y), 



1 

lef 



\ wej/ /z ' a / 



for a, 6 constant complex vector fields, and R satisfying 



\\R\\ 



LP 



< 



ICI 



\L\ \\(uJ 2 e fi I 8 + Q) k 3 



j,k=i 



li°°(flp) 



Furthermore, Y = (P — W l )Z is solution for (P + W)Y = in E and it reads 

Y=(0 ^H 1 | 7 1 / 2 E i )* 

with E, H solution for ([1]) in E. 

Proof: Let Y be defined by Y = (P — W l )Z, with Z the solution con- 
structed in Lemma [H If we denote 

we will prove that f 1 — f 2 — 0. 

Note that Y solves ffTTT) weakly, with f 3 solving the following decoupled 
equation: 

(-A + qj )f = 0, j = 1,2 

in the weak sense. 

Denoting L, R from Lemma [8] as 



L= ( I 1 (L 1 ) 1 | I 2 (L 2 Y )\ R=(r 1 (R 1 ) 



t I JZ 



(r 2 y y 



(25) 



the functions f 3 can be expressed as f ] = e^' x {rrij + Sj) with 

m x = C ■ L 2 - kI 1 si = --Dp ■ L 2 + (( + D - -Dp) ■ R 2 - nr\ 
m 2 = C ■ L 1 — kI 2 s 2 = ~Da ■ L 1 + (( + D — -Da) ■ R 1 - nr 2 . 

It is again a straight forward computation to check that 

(-A-2<- V + uj 2 e fi + qj)(m j + sj) = 0, j = 1, 2. (26) 
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Further, using (ITU]), PU|) and 



supp Da C -8(0; p), supp .D/3 C -8(0; p); 



one sees that Sj £ L^. Recall that equation f l2T)j) has a unique solution in L 2 ; 
whenever |£| > C(5,p) \\u 2 e fi + Qj\\ L oo( BD ), certainly it has to be the trivial 
one. Therefore, if irtj were also in . 
order to have rrij £ Lg , it is enough f 
be accomplished by choosing P,LP such that 



one. Therefore, if rrtj were also in L 2 , then p would vanish. Note that in 
order to have rrtj £ Lg, it is enough for its support to be compact. This could 



>. T 2 1/2 1/2,1 /• rl 1/2 1/2,2 

□ 

Next, we construct the same kind of solutions for the equation (P + 
Vy*)F = 0. Since (fT2l) holds and u 2 e po + Q has compact support in B{0; p), 
the same kind of arguments used in the proof of Lemma [8] can be carried out 
to state the following lemma. 

Lemma 10 Let 5 be a constant such that —1 < 5 < and let ( £ C 3 be 
such that C ' C = ^ 2 £oPo with 



lci>o-£ |(w%/i /8 + g) 
j,k=i 



L°°{B P ) 



for some constant C = C(5, p) > 0. Then, there exists a 



Z = e^ x (L + R) 

solution of (-AJ 8 + Q)Z = in E, with Z\ n £ H 2 (Vt; y), L — L(() constant 
and R = R(() satisfying 



R 



PR 



< 



C(S,p) 



~ ICI 
<o(ME 



(u 2 e poh + Q) k j 



j,k=i " 

8 

U 2 £ p I 8 + <2) 



j,k=l 









L°°(-B p ) V 







As a consequence of this lemma we get the following proposition, whose 
main difference with respect to Proposition [9] is the construction of solutions 
for the rescaled system instead of solutions for Maxwell's equations. 
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Proposition 11 Let ( G C 3 be such that C " C — ^^opo with 



\C\>C(p)J2\\(^oh + Q)* 

j,k=l 



L°°(B P ) 



Then, there exists a 

Y = e iC - x (M + S) 

solution for the equation (P + W*)Y = in E, with Y\ n e H l (Q; y), 

( C-a \ 



M 



ICI 



V Cx6 / 



for a, b constant complex vector fields, and S satisfying 



\ S \\L2(n-,y) < ^ 



■ J2 (||(w 2 £o/i /8 + Q)j 



r»(B p ) 



+ W. 



J IH,°°(fi) 



Proof: Let Z be the solution constructed in Lemma [TU| then by equation 
([TJ]), F = (P — W)Z is a solution of (P + = in E. Considering 

l = jjr ( o S* | o a* )* 

with a, 6 constant complex vector fields, the solution K can be expressed as 
Y = e^ x (M + S), with 

M = iP c L, S = PR + tP c R -WL- WR, 

where P^ is as in ffTTl) and S satisfies 



\\S\\ LH u;y) <C(6,p,n) J2 \\(u 2 eoHoh + Q) 

j,k=i 



(|i| + 


P 




L°°(-Bp) V 







+ c(5, P ,o)|L| ||(^ 2 ^o/8 + g)- 

j,k=i 



< 



L°°(n) 



C(O) 



L«(B P ) 
P 



ICI 



8 /ll 

^ (||(o; 2 eo//o/8 + Q)J 

.7,fc=l ^ 



+ w 



J IIL°°(Q) 



The last estimate is a consequence of Lemma [101 



□ 
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2.2 Proof of the log-type estimate 

The general ideas of this section go back to [I]. The most relevant difference 
is the use of a Carleman estimate. 

Let /ii,7i and //2,72 be two pairs of coefficients under the hypothesis of 
Theorem [2] and let us choose 

Ci = ~^ + i (r 2 + !^) V2 m + (r 2 + cA o/ , ) 1/2 rj 2 , (27) 

+ ^y / \i + (r 2 + c% / io) 1/2 r ?2 , (28) 



with r > 1 a free parameter controlling the size of |Ci| an d | C2 1 5 where £, 771, 772 



are constant vector fields satisfying \rji\ = 
for j = 1, 2. Note that Ci - C2 = -f and 



• ^1 , ^2 



C2 



1 772 1 = 1, r/i ■ r/ 2 = and 77^ • £ = 



Kil V2 \/2 



IC2 



y/2 V2 



Let us consider Z\ = e^ l ' x {L\ + Ri), Y\ the solutions stated in Proposition |9] 
corresponding to the pair /ii, 71 with > C(p, M). Recall that 



£1 



ICi 



/ Ci • oi \ 

1/2 1/2, 

Ci-&i 

\ 1/2 1/2 / 



l^il 



L 2 (n-y) 



< 



Kil ' 



Additionally, consider F 2 = e^ a ' E (M2 + 52) the solution stated in Proposition 
[TT1 corresponding to 1^2,12 with |£ 2 | > C(p,M). Also recall that 



Mo 



IC2 



/ C2 • a 2 \ 
-C2 x a 2 
C2 • h 
V C2X62 y 



C(p,fl,M) 
K2I 



Next we plug these solutions into the estimate ( TT8|) of Proposition [TJ with 
different choices of aj,bj. 

Choosing b\ = b 2 = and ai, a 2 such that 



1 y/2 V2 



■ CLi 



■ Vi r] 2 _ 
V2 V2} 2 ' 
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one gets, when r becomes large, that 

((Qi-Q2)z 1 \Y 2 ) n = 

= J e~^ x Qa(«! - a 2 ) + J(Va r Vai - Va 2 -Va 2 ) + (k 2 - dV 

+ 0((r 2 + |^| 2 )- 1 / 2 ), 

where the implicit constant in the symbol O is C = C(p,Q, M). Choosing 
o-i = cl% = and bi, b 2 such that 

i^= + ^=)-b 1 =(i^=+ V2 ^ r 



y/2 V2J \V2 y/2, 

one gets, when r becomes large, that 

((Q 1 -Q 2 )Z 1 \Y 2 ) Q = 

= J e-*" Qa(/3! - &) + i(V/3i-V/3! - V/3 2 -V/3 2 ) + (k 2 - k 2 )) ^ 

+ 0((r 2 + |e| 2 )- 1/2 ), 
where the implicit constant is C(p,Q, M). Write 

/ = In ~ a a) + ^(Vai-Vai - V« 2 -Va 2 ) + (« 2 - k 2 ; 

^A(ft - &) + ^(Vft-Vft - V/5 2 - V/3 2 ) + (/t 2 - k 2 ; 

where 1q is the indicator function of Q. By Proposition [7] and the properties 
of the special solutions, there exist three constants c = c(Q), C = C(p, Q, M) 
and C = C'(p, M) such that, for any r > C one has 

1/(01 + l?(0l < C7 (BCfcCCx^a))^^ + (r 2 + KI 2 )- 1 / 2 ) . 
Note that, for s± < and R > 1, one has 

k ( E) + ibiik(E) = / (i + io 2 r (i/cor + mi 2 ) d$ 

At\<R 

+ / (i + iei 2 r(i/(oi 2 + i?(oi 2 K 

<C( J B(^(C 1 ,C 2 ))e c ^+ r- 1 ) 2 Al + |r| 2 ) s V 2 dr 

Jo 

+(i+R 2 r(\\f\\l m + 
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Therefore, 



H s i(E) 



+ \\9\\ H s 1{E) < C (B(5 C (C U C 2 ))e« R+ ^ + r-'R 3 ^ + R s >) . 



Now we choose R in such a way that r 1 R 3 / 2+s ^ decays as R Sl , that is, 
R = r 2//3 , hence 

11/11^(0) + IMIn^n) ^ C ( s (^(Ci, C 2 ))e CT + r 2 / 3 -) . 
On the other hand, the a priori bound was chosen to have 

\\f\\H»Hn) + \\9\\Hs 2 (n)< c ( M )i 
for < S2 < 1/2. Finally, by the interpolation estimate 



there exist two 

constants C = C'(p,M) and C = C(p,Q, M,u) such that, for any r > C, 
the following estimate holds 

\\f\\»V) + WdWmn) < C {B(S {C 1} C 2 ))e~ + r^) e , (29) 

with = 9s 1 + (1 - 0)s 2 . 

The idea now is to transfer this estimate from /, g to the difference of 
the coefficients n\ — n% and 71 — 72. This can be accomplished by using the 
following Carleman estimate. 

There exists a positive constant C(Q) such that, for all h < 1 and any 
function G C 1,1 (fi), the following estimate holds 



L 2 (H) 



3¥ >//i 



V0 



L 2 (n ; C 3 ) 



< 



< C (h^\\e' p/h A(j)\\ 2 L2{n) + h\\e ip/h 



\L 2 (an) 



V0 



L 2 (9n ; C 3 ) 



where </> = 1/2| X — Xq 1 2 with Xq £ f2. The constant here depends on the 
distance from Xq to Q and on the diameter of Q. A Carleman estimate of 
this type can be found in [TO]. 
A simple computation give: 



/ = ln7i 



9 = In/^i 



1/2 



1/2 



a / 1/2 1/2 \ . / 1/2 

A (7i -7 2 ) +9/(71 

A / 1/2 1/2% . / 1/2 

-lh ) + Wi 



7 2 ) +P/W -^2 
^2 /2 ) +P 9 (7i V2 -7 2 1/2 ) 



1/2 



where 
9/ = - 

% = - 



A 72 /2 , 2 1/2/ 1/2 . 1/2 

+ w 71 (7i /^i + 7 2 ^2) 



^ 1/2/ 1/2 1/ 

P/ = 7i7 2 Oi + ^2 

2 1/2/ 1/2 l/2x 

p 5 = -u nin 2 ' (7/ +72 )■ 



72 

A/i 



1/2 



1/2 



. 2 1/2/ 1/2 1/2 \ 

w //a' (//i 71 + fa 72) 



1/2 
M 2 
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Note that, thanks to the a priori bounds, we have the following differential 
inequalities: 

|A( 7l 1/2 - 7 2 V2 )I < C(M)(\f\ + | 7l 1/2 - 72 1/2 | + \^' 2 - /if |), 
|A(/if - /if )| < C(M)(\g\ + | 7l 1/2 - 72 1/2 | + |^ /2 - /if |). 

In order to simplify the notation, let us write 0i = 7 f — 7 f and 2 = 
/if — /if. By the differential inequalities written above and the Carleman 
estimate, one has 



< C" £ P /fc fc||L ( n) + * ll e ^llL(«,) + ^ 11^0,11^^3)) 
+^ 4 (||e" /fc /||L( n ) + ll^||L(n)). 



i=i,2 



where the constant is C" = C"(Q,M). The terms /i 4 || e^/ 71 ^ | |^ 2 m)' w ith 
j = 1, 2, can be absorbed by the left hand side of the inequality. Hence, if 
di = mi{d e (x] x ) 2 : x G Vt\ and d 2 = sup{<i e (a;; x ) 2 : x G f2} we get, for any 
/i < C'^M)" 1 / 3 , that 



e dl/h E ( fc II^H^W + ^ ll V ^ll^(^;C3)) < CV 

J"=l,2 

^ 4 (ll/ll^n) + IMk (n) ) + E ( h + ^ H V< ^I W 3 )) 

J"=l,2 

But now we can easily estimate 

\\M LHd n) < Chi -7ilL-(«i) < CS(fc(Ci,C 2 )), 
l|V0i|| L2(an . C 3 ) <c(|| 7 i + ||V( 7 i-72)|| L 

<CS(5 c (Ci,C 2 )), 

II0 2 || LW < cii^ -HU^ < c^Md,^)), 

llV^H^^.cs) < C (||/ii - ^Hloo^) + ||V(/ii - Ai2)|| L oo (a n;C3)) 
<CB(6c(C lt C 3 )) 1 
hi - 72|| L2(n) + ||V(7i - 72)11^2(^3) < C M|0i|| + l|V0i|| LW J, 

ll^i - V2\\ L 2(a) + ||V(/ii - A*2) Ilx,2(n ; c3) ^ C (ll02|| L 2( Q ) 
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The constants above depend on the a priori bounds M. These inequalities 
and estimate f[2"9j) gives us 

Il7i - TiH^oj + Ik - ^11^(0) <Ce^ (B(8 C {C X , C 2 ))e CT + r 2 ^)^ 

+ Ce^B{5 c (C l ,C 2 )), 

where d 2 > d x , si < < s 2 < 1/2, c = c(fi), C = C(p, tt,M),r> C'(p, M) 
and h < C"(Cl, M) -1 / 3 . To end up with the estimate given in the statement, 
it is enough to note that 

2 sxs 9 2 

6 S 2 — Si 3 

and to choose the parameter r as 

r=-ilogS( < J c ((7 1 ,C 2 )). 

A Appendix 

In this appendix, we state some elementary concepts and collect some known 
facts relative to Lipschitz domains and Sobolev-Besov spaces. 

A.l Lipschitz domain 

Definition 12 Let Q be a nonempty, proper open subset of E and fix a point 
Xq on its boundary dfl. We say that Q is a Lipschitz domain near xq if there 
exist 

(i) a plane q C E passing through x an d a choice of a unit vector iV 9 
normal to q; 

(ii) some euclidean coordinates £ : E — > IR 3 such that £0(^0 ) J = for 
j = 1,2,3 and 8 (x) G M 2 x {0}, for any x E q (for short, we shall 
denote £o(a ; )' : ' by y- 7 ); 

(hi) and an open cylinder C%° = {x G E : \y'\ < c%, \y 3 \ < c 2 } -called 
coordinate cylinder near x$- such that 

cz C2 nn = c*° C2 n {1 e e : j/ 3 > y 2 )}, 

C*° C2 ndQ = C'Z C2 n{xEE:y 3 = <f)(y\ y 2 )}, 
C^ C2 n H c = C*° C2 n {* G E : y s < <P(y\ y 2 )}; 
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for some Lipschitz function <p : M. 2 — > R satisfying 

0(0) = 0, and |0(2/V)|<c 2 if \y'\ < c x . 

Finally, we say that Q is a Lipschitz domain if it is a Lipschitz domain 
near every point x G <9fi. 

In this definition the superscript c denotes the complement of a set, relative 
to E; and \y'\ 2 = \y Y \ 2 + \y 2 \ 2 . 

Recall that for a Lipschitz domain there exists a measurable unit normal 
vector field N along dQ pointing outward. In this context, we can set the 
following integration by parts formulas 

- [(V-u)fdV= [ u-VfdV- I (N-u)f\ dn dA (30) 
Jn Jn Jan 



and 



I (Vxu)-vdV = / u-{Vxv)dV+ / (N xu)-v\ dn dA, (31) 
Jn Jn Jan 



where / is a smooth function on Q and u, v are vector fields on Q. Here dA 
stands for the area element and .\ga denotes the restriction to dfl. 

A. 2 Sobolev and Besov spaces 

Most of the facts collected here can be found in 1 1 1 1 and the references therein. 



Definition 13 For any s 6 R, define the potential Sobolev space on E as 

H S (E) = C§°(E) IMIhS(e) , 

with the norm 

imiW)= / (i + i£i 2 ri/(oi 2 ^. 

When s = this space is L 2 (E). 

Denote by f\ n the restriction of / G L 2 (E) to H C E. 
Definition 14 For any s > 0, define the potential Sobolev space on Q as 

H S (Q) = {f\ u :fe H S (E)}, 

with the norm 

IbiU(n) = inf {ll/ll^(E) : /In = g}- 
28 



On the other hand, for sGK, define 

H°(Q) = {/ G H S (E) : supp / CO}, 

with norm 

Il/lli?g(f2) = H/H# s (E) • 
Finally, for s > 0, define the space H~ s (fl) as the dual of Hq(Q), that is, 

It is well known that C°°(n) = {f^ : / G C°°(E)} is dense in H S (Q), for 
sei When s = this is L 2 (fi). 

On the other hand, H S (K) is a complex interpolation scale for s G R; that 
is, for any s\, S2 G R, one has 

[H Sl (E),H S2 {E)} 6 = H S (E), 

with s = 9s\ + (1 — 6)s2 with 9 G (0, 1). Moreover, the extension by zero 
outside Q is a bounded linear operator from H S (Q) to .£P(E) for —1/2 < s < 
1/2 (see [22]), which allows to identify the spaces H S (Q) and Hq(Q,) for the 
same range of s. With these facts in mind, notice that 

ll/lk^)<^ll/ll^(n)ll/ll^ ( n)> ( 32 ) 

for s = Bs x + (1 - 6>)s 2 with s x G R and -1/2 < s 2 < 1/2. 
Recall that 

W 2 '°°{Q) = {/ G : 9°/ G < \a\ < 2}, 

with the norm 

ll/llw 2 '°°(n) = H^/IIl 00 ^) • 

0<M<2 

We next define the Besov spaces in an intrinsic way. To do it we introduce 
the functional 

defined for / G 5(R 2 ) -the space of rapidly decreasing functions. 
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Definition 15 For < s < 1, let us define the Besov space 
B S (R 2 ) = {f e L 2 (R 2 ) : I s (f) < +oo}, 

with the norm 

ll/W)HI/IL w + n/)- 

Now we shall extend these Besov spaces on M 2 to Besov spaces on dQ. Let 
xi,...,x n belong to dQ and Yi, . . . , Y n be T 3 ; — Cc?, C2 H dQ for j — 1, . . . , n, 
such that <9fi = Ti U ■ ■ • U T n ; and consider a partition of unity X\i - ■ ■ iXn 
subordinate to r 1; . . . , T n . We shall say that / G B s (dQ) for < s < 1 if 

( Xj f)oEr%<l> j (.))eB'(R 2 ), 

for any possible choice of points and any partition of unity related to them 
as above. Here £j and <f>j are, respectively, the euclidean coordinates and the 
function defining the boundary locally, corresponding to the point Xj. The 
norm defined on these spaces will be given by 

n 
3=1 

Xj G dQ, supp(xj) C Fj, j = 1, . . . , n}. 

One of the reasons to introduce these spaces is to describe the properties of 
the trace operator. Namely, the trace operator .\qq : H s (£l) — > 5 s_1 / 2 (<9r2) 
is well-defined, bounded and onto, whenever 1/2 < s < 3/2. In addition, it 
has a bounded right inverse whose norm is controlled by s and the Lipschitz 
character of f2. 

For later references, we give the following lemma. 

Lemma 16 Let s, e be such that < s < 1 and < e < 1 — s. Then, 
there exists a constant C(s, e) > such that, for any g G C°' s+t (d£l) and any 

feB s (dn), 

Remark: The constant C(s, e) given here blows up when e becomes 
small. 

Proof: Let x±, . . . , x n belong to dQ and . . . , r„ be such that Tj = 
Cci lCa n dfl for j = 1, . . . , n, such that dQ = I\ U • • • U T n . Let A > be the 
Lebesgue number associated to {r.,}™ =1 and define 

f, = {i 6 T, : inf d e (x, x') > A/4}. 
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Note that dQ — f x U • • • U f n . Let us consider a partition of unity Xu ■ ■ ■ > Xr. 
subordinated to T 1 , . . . , T n and denote 



Here lr 3 stands for the indicator function of Tj. Consider Xk G Tj with 
k = 1,2, we can write x^ = E~ x {y^ 4>j(yk)) for G M 2 . By the Lipschitz 
character of <pj, one has 

d e (x 1 ,x 2 ) < C\yi - y 2 \. 

Then, noting that 

(I s (a f )f < 2 sud \9j(y' + y)-9j(y')\ 2 n f |,2 /" 1 ... 

V \93J3)) SUp L,|2( S +e) IIJJIIL2(R2) / U,|2(l-e) a f 

2/'6£ 3 (r 3 ) IS/1 v ; Jc|y|<A/5 12/1 v ; 

C[2/|<A/5 

+8||^|ii 00(Ra) ||/iii^ (R3) / l, Ji+s) rf y + 2 II^-IIx~cr») ( /s (/?)) 2 > 

^C|y|>A/5 

we can achieve the result. □ 

Definition 17 For < s < 1, define the space B~ s (dQ) as the dual of 
B s (dQ), that is, 

£- s (<9ft) = (B s (dn))*. 

In the same conditions as in Lemma [161 one nas > by duality -recall that 
( g f\h) = (f\gh) for any h G 5 s (<9ft)-, that 

ll#/II.B-*(an) — C \\9\\c > s +c(dn.) Wf\\B- s (dn) • 

The constant here is the same as in Lemma [TBI and blows up when e becomes 
small. 

Definition 18 For s G R define 



where 

3 

11 II 2 — II 0')ll 2 

\\ U \\Hs(n;C3) - W U \\H»(Cl)- 

i=i 

Here AfE|jy = {u\q : u G XK}. When s = this space will be denoted by 
L 2 (ft;C 3 ). 
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Definition 19 For < \s\ < 1 define 

B s (dn-,c 3 ) = xE\^ UBS{dn -- c3) 

where 

3 

II II 2 — II (i)ll 2 

W w \\B°(dn;C3) - \\ w \\B"(dn) 

i=i 

and AfEl^n = '■ u G XK}. In the same way, let us define 



L\dn-X A ) = *E| an " 

where 

IMli^(an;c») - H w "' IIl 2 (9Q) 



2 = ^IU„0')ll 2 

i=i 



Let us point out that, for < \s\ < 1 and < e < 1 — \s\, 

\\9 w \\B»(dn;C 3 ) — C lbllc°>n+ e (9n) IMIb^scjC 3 ) ■ 

The constant above is the same as the one in Lemma [16] and, once again, it 
blows up when e becomes small. 

On the other hand, by using the trace operator component by component 

.\ dn : H S (Q; C 3 ) — > B*' 1 ' 2 ^; C 3 ) 

is bounded and onto, whenever 1/2 < s < 3/2. 
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